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The complex amplitudes of the electronic wavefunctions on different sites are used as 
Kramers variables for describing Electron Transfer. The strong coupling of the electronic 
charge to the many nuclei, ions, dipoles, etc, of the environment, is modeled as a thermal 
bath better considered classically. After elimination of the bath variables, the electron 
dynamics is described by a discrete nonlinear Schrodinger equation with norm preserving 
dissipative terms and Langevin random noises (at finite temperature). 

The standard Marcus results are recovered far from the inversion point, where atomic 
thermal fluctuations adiabatically induce the electron transfer. Close to the inversion 
point, in the non-adiabatic regime, electron transfer may become ultrafast (and selec- 
tive) at low temperature essentially because of the nonlinearities, when these are ap- 
propriately tuned. We demonstrate and illustrate numerically that a weak coupling of 
the donor site with an extra appropriately tuned (catalytic) site, can trigger an ultra- 
fast electron transfer to the acceptor site at zero degree Kelvin, while in the absence of 
this catalytic site no transfer would occur at all (the new concept of Targeted Transfer 
initially developed for discrete breathers is applied to polarons in our theory). 

Among other applications, this theory should be relevant for describing the ultrafast 
electron transfer observed in the photosynthetic reaction centers of living cells. 



1. Introduction 

According to transition state theory, chemical reactions decompose into elementary 
reactions among which electron transfer is ubiquitous 0. The time required for an 
electron transfer (ET) between different molecules is expected to be minimum at 
the so-called Marcus inversion point in the space of parameters. However, it is 
precisely in the vicinity of this inversion point that the standard adiabatic (Born- 
Oppenheimer) approximation used in Marcus theory, breaks down. We propose a 
non adiabatic dimer model which should improve the Marcus theory in this regime. 
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We also show that our approach extended to trimer models predicts the possibility 
of catalytic ultrafast electron transfer when appropriate conditions are fulfilled. 

The Marcus theory □ well describes in many cases the features observed for 
ET between two molecules. It is essentially an adiabatic theory where the atomic 
fluctuations are supposed to be slow at the scale of the characteristic time of the 
electron dynamics so that the wavefunction of the electron may remain practically 
an eigenstate of the time dependent potential created by the atoms. This adiabatic 
assumption is valid when the largest phonon energy hu> p h remains much smaller 
than the smallest excitation energy of the electron. In a two level system, this 
energy is essentially the distance Tiuj e i between the two levels. This condition is 
very often fulfilled in real systems because electronic energies are generally much 
larger than phonon energies. 




Fig. 1. Free energy versus Reaction Coordinates of the system Donor- Acceptor when the electron 
is on the Donor (top left curve D) or on the Acceptor for several redox potentials in the normal 
regime (top right curve), at the inversion point (middle right curve) and in the inverted regime 
(bottom right curve). The chemical reaction energy is the distance between the energy minima 
AG . The energy barrier is AG*. The electronic excitation energy on the Donor at fixed Reaction 
Coordinates is A e ;. 



Since in the adiabatic case the electron dynamics is essentially driven by the 
dynamics of the atoms, ET is essentially induced by their fluctuations. The free 
energy of the system as a function of the reaction coordinates (describing the global 
atomic configuration), depends on the electronic state which could be either on 
the Donor molecule or the Acceptor molecule. The well-known Marcus scheme is 
represented in fig.[l] (see ref.i for details). In the normal regime, the ET from the 
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Donor to the Acceptor at fixed Reaction Coordinates requires a positive energy 
hu> e i = — A e i. ET is a thermally activated process since thermal fluctuations of the 
lattice are necessary to overcome the energy barrier AG*. In the inverted regime, 
this energy — A e ; is negative. Although ET could be achieved at low temperature 
by a photon emission at frequency (energy) fuv e i = A e ; (photoluminescent chemical 
reaction 13), activation processes above the energy barrier AG* become by far more 
efficient and prevalent at higher temperature. 

AG* turns out to be just zero at the inversion point when A e / = (see fig.0). 
This is the regime where ET is expected to be at maximum speed and still effective 
at low temperature because of the absence of energy barrier. However, in the vicinity 
of this inversion point the validity of the adiabatic hypothesis necessarily breaks 
down, since the characteristic energy |A e ;| of the electronic excitation becomes 
small. Our approach improves the theory for ET to be valid in this non-adiabatic 
regime as well. 

2. A Non-adiabatic Model for ET 

We consider a single electron tunneling between Donor (D) and Acceptor (A) sys- 
tems representing large molecules with many vibrational degrees of freedom (see 
fig.|j). Each of these molecules a (a = D or A but more molecules may be involved) 
is supposed to involve for simplicity a single electronic state with a wavefunction 
\^f a >= *& a (r; {uf }), where r is the space coordinate, uf the phonon coordinates, 
and we assume that the Born-Oppcnhcimcr (adiabatic) approximation. This ap- 
proximation is valid when the other electronic states on this molecule are far apart 
in energy from the considered state \*& a > at the scale of the maximum phonon 
energy. Then this electronic wavefunction can be considered as a function of the 
molecule phonon coordinates {uf}, including possibly those of the environment and 
in particular, the solvent. 

Within a standard tight-binding representation, the state of the electron in the 
whole system, has the form X^V'al^a >■ We use as Kramers reaction coordi- 
nates the complex amplitudes ip a of the electronic wavefunction. We have the 
normalization condition ^ a |i/'q| 2 = 1- 

It is convenient to define first the minimum energy of the system of the two 
coupled electronic states HT({ip a }) at fixed collective variables {ip a }- Then, the 
Hamiltonian of the interacting electron-phonon system can be written as 

h = iMOM) + £fl£ h ({hM 2 ,«i a ,p?}) (i) 

a 

The electronic density \tp a \ 2 on each molecule a couples to the coordinates of the 
same molecule assumed to be harmonic and thus consisting of a collection of in- 
dependent harmonic oscillators i with position-momentum coordinates u Ql i,p Q .i, 
mass m a ,i, and frequencies io a ,i 

H* h ({\il) a \ 2 ,u a ,i,p a ,i}) = (u a .i - k a .i\tp a \ 2 ) 2 + 1 7——p%j (2) 
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Electronic 




Phonon bath 



Fig. 2. Schematic representation of a Donor- Acceptor pair with electronic levels interacting with 
phonon baths in the normal Marcus regime 

In principle this coupling energy involves all possible interactions with the atomic 
coordinates and in particular, the chemical energies and the electrostatic energies. 
The latter ones could be especially important in biomolecules which are polyelec- 
trolytes surrounded by ions and highly polarizable water. 

In general HtH^o}) has not the Hermitian form < {ip a }\HT\{ip a } > where 
Ht is a linear operator but is highly nonlinear as a consequence of the electric 
field and the molecule reorganization generated by the density variations of ip a . It 
is convenient to split this Hamiltonian in several parts 

#t({iU) = ]T H a (\^ a \ 2 ) + ff/({K| 2 }) + H t ({^ a }) (3) 

a 

where H a (\ijj a \ 2 ) is the energy of the isolated molecule a which depends only on its 
electron density I a — \ip a \ 2 ■ It is sufficient to expand this energy at second order 
in electronic density for obtaining the main physical features 

H a (\^ a \ 2 )=^ a \ 2 + lx a \^ a \ 4 (4) 

\i a is the linear electronic level at zero occupation. \a — Xa + Xa ^ s the sum of 
two contributions. Xa * s the positive coefficient for the energy of the electric field 
generated by the charge I a without lattice reorganization (capacitive energy). This 
coefficient takes into account the electronic dielectric constant e^. 

Xa lS the negative coefficient of the energy gain from the local reorganization 
due to the presence of the electron. This energy involves for a part chemical bond 
energies which could be broken or created and electrostatic terms involved in the 
static dielectric constant. Actually, in our model (||) the reorganization energy can 
be explicitly calculated and we get Xa — ~ J2i m a,i LO a i^a %■ 
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It is essential to remark that that the sum of the two contributions \a might 
be positive or negative depending whether it is the electrostatic energies or the 
chemical energies which are prevalent. We expect for example that when the elec- 
tronic state a belong to the inner shell of a transition metal ion (which could be 
embedded in a large biomolecule), the electrostatic energy is prevalent so that Xa 
is positive. When it belongs to a chemical bond (or a ring of bonds), it is more 
likely negative. Actually, only ab initio calculations could estimate the real values 
of these coefficients. The essential physical consequence of the nonlinearity is that 
the electronic level on molecule a 

E a = = [l a + XalV'al 2 (5) 

Ul a 

depends on its electronic occupation density. 

Hf({\ip a \ 2 }) in is a small extra term due to molecule interactions which 
which could be for example the Coulomb interaction energy Hf({\i/j a \ 2 }) = 
J2 a a Ca,p\ipa\ 2 \ipf3\ 2 where C Qj /3 are mutual capacitance coefficients. 

H t ({ip a }) in (|J) is also a small extra energy term due to orbital overlaps but it is 
essential because it allows ET. We may choose for simplicity the form H t ({ip a }) — 
\ a pi/ia^p + c.c. The transfer integrals X a ,p & r e assumed to be small (compared to 
differences of electronic energy levels) and may be of the order of phonon energies. 

In an ideal anti-adiabatic regime (but non-realistic) where the transfer integrals 
X a ,(3 between different molecules would be much smaller than all phonon energies 
hu) a .i, ET would be much slower than the phonon dynamics. Then, the phonons 
could be eliminated as fast variables following adiabatically the slow electron vari- 
ables and Hrdipa}) would be the exact Hamiltonian describing the electron dy- 
namics through the set of nonlinear Hamilton equations iKip a = OHt / 'dip a . 

Actually, nonlinearities may generate energy barriers and even in the absence 
of energy barriers, the absence of energy dissipation does not allow ET to a lower 
energy level. Ultrafast ET requires an efficient energy dissipation which can be 
obtained only by interaction of the electron with a phonon bath in the regime 
intermediate between adiabatic and anti adiabatic ^. 

Usually, the lattice reorganization due to the presence of an electron on a 
molecule is large and involves the coherent creation of many phonons. It is thus 
legitimate to treat this phonon bath classical while the electron dynamics remains 
quantum. The dynamical equations of the coupled system (jll@) are 



ihip a = — — - } m a ik ati uj a i(u a ,i - k a ^\ip a \ )ip a (6) 

U a .i + wiA u a,i ~ k a ^\lp a \ 2 ) = (7) 



a The dynamical coupling of the electron with the electromagnetic field also generates energy 
dissipation. However, this coupling which is weak has to be treated quantum using Fermi golden 
rule which yields relatively long life time to quantum excitations. We neglect here this energy 
dissipation since we are interested in ultrafast ET. 
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The harmonic motions u a ^(t) can be explicitly obtained from the linear equa- 
tions (Fij) □ as the sum of functions of the time dependent driving force |7/i Q (i)| 2 
and a solution of the equation without driving force. Actually, this term physically 
corresponds to thermal fluctuations of u a ^ and thus is random. Then, substituting 
Ua,i{t) hi eq.(§) yields the fundamental equation for non-adiabatic electron dynam- 
ics (which preserves the norm \ipa\ 2 ) 

ih ^ a =^+fj t T a (t-T)^-dT+ut))^ (8) 

where T a (t) = J^. m Q ^w 2 Jt^ i cos (w a ^t). If there are many phonon modes with a 
rather uniform distribution, T a (t) can be assumed to be a smooth decaying function 
of time. It generates energy dissipation as a kernel in eq.(||) (the absorption rate in 
energy of a charge fluctuation at site a at frequency lo is nothing but the product 
of the square of its amplitude with the Fourier transform of T a (t)). We also have 
r a (0) = — Xa- The time dependent potential Ca{t) 1S produced by the thermal 
fluctuations of the lattice. It is a colored random Langevin force with correlation 
function which fulfills < ( a (t + t)C q (t) > r = T a (t) ksT at temperature T a. 

Thus, the effect of non-adiabaticity is to transform the standard linear Schrodinger 
equation describing the dynamics of the electron into a nonlinear Schrodinger 
equation (||) with norm preserving energy dissipation terms and with random 
colored time dependent potentials generated by atomic thermal fluctuations. 

3. ET in the Dimer Model 

We first show that far from the Marcus inversion point we essentially recover the 
basic result of the standard theory □. The initial Hamiltonian (0J^) restricted to a 
dimer model (a = D or ^4) readily yields the energy surfaces schematically shown in 
fig.|l|. Neglecting the small interaction energy terms between the molecules we obtain 
the essential parameters of this theory which are AG = [in + \xd — — \xa, 
A el = AG + i( x g + Xa 1 ) = Mb + xg - Ma and AG * = -^li/Wxl + X^)) whcrc 

Ma = Ma + Xa A 

The energy variation Et{Ia) = Hn(l — Ia) + Ha(Ia) — Hd(1) of our dimer as 
a function of the electron density I a = \iPa\' 2 on the acceptor is 

E T {I A ) = - Mb - Xd)Ia + ~(xd + Xa)I\ (9) 

and — .E't(I) = AG is the chemical reaction energy. There is always an energy bar- 
rier between Donor and Acceptor when ^p(O) > or equivalently (ia < M-D + Xd- 
Otherwise, the minimum of energy is not necessarily obtained for a total transfer at 
the Acceptor when xd + Xa > 0. The derivative dExIdlA = Ea — Ed (see eq.((|)) is 
the difference of the electronic levels on the Acceptor and the Donor at the transfer 
I A- Thus, resonance between Donor and Acceptor implies a zero derivative. 

For recovering the same results as the Marcus theory from our equation, it is 
essential to note that the phonon spectrum has a cut-off at relatively small frequen- 
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cies u) c . Beyond this frequency, the Fourier spectrum of T{t) is zero. Thus, when 
the characteristic energy of the electron dynamics, which is the energy difference 
Ed — Ea between the electronic levels, eo.(||), is larger than the phonon energy 
huj c , there is no more energy dissipation □. Then ET cannot be achieved at zero 
degree K but requires thermal fluctuations. Actually, this is the regime of validity 
of the adiabatic approximation. 

The electron density which is initially on the Donor \4>d\ 2 = 1 and \ipA\ 2 = 
remains practically constant. Eq.@ gives in this case the random potential as due 
to phonon variables, Ca(t) = — J2i m a,ik a ,i^a i( u a,i — k a ^\ip a \ ). The slowly varying 
potential Cd(^) makes the energy level of the electron E D (t) = ^d + Xd+ Co(t) on 
the Donor time dependent. The unoccupied energy level on the acceptor £U(£) = 
fiA + Ca(£) fluctuates similarly. It occurs statistically that ED(t) « EA(t) or 

+ XD ~ »A + CdW - Ca(*) » (10) 

induces an almost resonance between Donor and Acceptor so that the electron 
could tunnel (see figJJ). If we discard the details concerning the probability of this 
tunneling process and neglect its intrinsic time which is generally short compared 
to the characteristic time for reaching the resonance, the characteristic time for ET 
is mostly related to the time required to reach the resonance []. 

Condition |l0| may be compared with the condition for the intersection of the 
two free energy surfaces shown in fig.|l|, which can be written as 

^'d+Xd- + Cd(*) - U(t) = A e , + ( D (t) - U(t) = (11) 

This condition is similar but different from our resonance condition (^) . The reason 
is that the energy level variations due to the Coulomb energies are not taken into 



b The electron tunneling problem has been considered differently in the literature which disti: 
guishes between adiabatic processes (strong reactants) and diabatic processes (weak reactants) 
Our non-adiabatic theory (potentially) describes both cases as well as the intermediate cases. 
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Fig. 4. Sketch of the electronic level fluctuations in the normal regime (left) and in the Marcus 
inverted regime(right). 

account in the Marcus theory unlike the reorganization energy ( if Xa = conditions 
( [To| ) and (|ll]) become identical). Nevertheless, we can also interpret the probability 
of reaching the resonance as shown in figj^a in terms of an activation process with 
an activation energy A'G* = -A' e1 /(2(x 1 d + Xa)) AG * where A' el = A el + (xg + 
Xa)/2 = A*£> + Xd — A< a is different from A e i because of the Coulomb terms. 




time (ps) time (ps) 

Fig. 5. Electron density on the Donor and the Acceptor versus time for the dimer model at the 
inversion point and zero degree K where (ip = 2, \D = — b MA = b XA = —0.75, \ad = 10 -2 , 
7o = = 1. (left) or 20. (right) (the time unit is 1 ps=10~ 12 s for energy in units of eV) 

Nevertheless, our approach confirms the existence of an inversion point when 
resonance is obtained at zero degree K for (d = Ca = 0. Then, (|l0|) yields A' el = 
or fiD + xd — (J-A which again is different from the condition A e ; = in the original 
Marcus theory. 

At our inversion point, the energy profile Et(Ia) has a zero derivative at the 
origin I a = 0. Fig. [5] shows that there is no energy barrier only when \d + xa < 
R. We check that the initial resonance Ed = Ea triggers ET. If the derivative 



c The only case implicitly considered in the Marcus theory is for negative \ a = x„ < 0. 
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dEx I dlA remains smaller than htu c which is equivalent to a small reaction energy 
AG < h,uj c /2, ET can be achieved at zero degree K without thermal fluctuations. 
Fig.|| shows two examples. The speed of ET strongly depends on the damping 
which is related to the coupling to the phonon bath. In these examples and in the 
following we chose for simplicity Hf({\ijj a \ 2 }) = in (||). If the phonon frequency 
cutoff oj c is large compared to the characteristic electronic frequencies, a reasonable 
approximation is to assume that T(t) = 2j a 6(t) is a Dirac function. Then, 
ftoo r(£ — T ) d ^ dr w 7q {t) in eq.(||). There is an optimal damping constant 
(7a ~ 40) where the characteristic time required for ET is minimum. 

ET is triggered at zero degree (but slows down) when escaping only on one 
side of the inversion point when fiA < /J-d + Xd- On the other side, it is blocked 
at zero degree K because of the appearance of an energy barrier. However it is 
complete only when dEx / dlA remains always negative with a modulus which never 
exceeds the phonon cutoff energy %ui c (for having efficient energy dissipation). In 
summary, our approach yields results which are qualitatively similar to those of the 
Marcus theory far from the inversion point but with a redefinition of characteristic 
parameters. It yields more detailed features not predicted by the original Marcus 
theory, close to the inversion point. 

4. Principle of Catalytic ET in a trimer model 

ET could be fast for a Donor- Acceptor system only in special conditions close to the 
Marcus inversion point and when the chemical reaction energy is small compared 
to the phonon energy cutoff. 

We now show that we can take advantage of a third catalytic site weakly coupled 
to the Donor for triggering at zero temperature an ultrafast ET from the Donor to 
the Acceptor while in the absence of catalyst, a large energy barrier would prevent 
any transfer at zero degrees. 

It is clear that ET could become fast only in case of resonance or almost res- 
onance within the phonon energy range. Otherwise, the electronic level E a on a 
molecule depends on its occupation density E a = dH a {I a ) / dl a = fj, a + Xa T a 

How to get complete ultrafast ET between Donor and Acceptor at zero degree 
K not at the Marcus inversion point? We suggest to take advantage of the nonlin- 
earities for inducing relatively slow oscillations of the electronic level on the Donor 
(or on the acceptor) which could produce a resonance between Donor and Acceptor. 

The simplest way is to obtain these energy oscillations by thermal fluctuations 
as shown figj^. This is the standard situation described above where Marcus theory 
is recovered. 

Another way is to induce artificially these electronic level oscillations by exciting 
a specific phonon well coupled to the electronic level. Note that this is what happens 
systematically when creating an exciton by the absorption of a photon before the 
Franck-Condon relaxation. This situation very likely occurs in some real systems 
but we shall not discuss it here. 
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Donor 




Catalyst 

Fig. 6. Sketch of electronic level oscillations on the donor system inducing resonance with the 
acceptor level. 

For inducing these level oscillations, we propose to use a third site, a "Cat- 
alyst" , which is appropriately tuned on the Donor. For that purpose n we use the 
phenomenon of Targeted Energy Transfer (or Targeted ET in this case) ErQ. A Donor 
and an Acceptor at the inversion point are in resonance but in general they do not 
remain in resonance during the transfer. However, there is a specially interesting 
case, when \d + XA = (which is exactly solvable i) for which the resonance 
persists all along the transfer. Since cIEt{Ia) / 'cIIa = 0, Et{Ia) = 0, the reaction 
energy is zero. As shown in I in the absence of damping, the electron slowly os- 
cillates between the Donor and the Acceptor, which is now the Catalyst (C), with 
a frequency corresponding to the transfer integral Xcd- The electronic level oscil- 
lates with the half period. This situation requires to associate a soft electronic level 
with an appropriate and well defined hard electronic level ( which could involve a 
metallic ion as suggested above) (fic = Hd + Xd and xc — ~Xd)- When there is en- 
ergy dissipation, this oscillation is damped and converges to the covalent state with 
equal density on Donor and Catalyst (see fig.0) . Indeed, the small transfer integral 
raises the degeneracy at zero coupling. The range of variation of the energy level, 
which is the interval [fijj + xdjI^d] in the undamped case, is reduced to the inter- 
val [/i£> + xi)/2,/Lt£)] in the overdamped case. However, this binding energy is very 
weak and negligible since Xcd is small. Then, small thermal fluctuations become 
sufficient to generate giant charge fluctuations. In practice, Donor and Catalyst do 
not bind chemically but they could trigger ultrafast ET to another molecule. 

We now test the principle of catalysis suggested by figj|. We consider a trimer 
model with Hamiltonian 

H tr = HdHd\ 2 + ^XdI^dI 4 + MclV'cl 2 + ^Xc#c| 4 + ^a\^a\ 2 + \xa\^Pa\ A 

+Xcd{^ d ' 1 Pc + cc) + Xca^a^c + cc) + Xad^d^a + cc) 

We assume that Donor and Acceptor are both soft (xd < 0,Xa < 0). The transfer 
energy is positive which means ha + Xa/2 < hd + Xd/2- 
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In the absence of Catalyst (Xcd = ^ca — 0) there is a large Kramers energy 
barrier between Donor and Acceptor, which implies /j,d + Xd < Ha- Thus, we are 
in the normal Marcus regime. ET is impossible at zero degree K and remains slow 
up to relatively high temperature. 

We now introduce the Catalyst. Donor and Catalyst are tuned for Targeted 
Electron Transfer which requires xc = ~Xd and He = Hd + Xd- The initial 
electronic level ha of the acceptor should belong to the variation interval of the 
electronic level of the Donor-Catalyst system which yields hd + Xd < HA < Hd 
at weak damping or hd + xd/% < Ha < Hd at strong damping. Fig.|| sketches 
Hd(Id), Hc(Ic) and Ha(Ia) when all these conditions are fulfilled. 

Fig.|| shows that under these conditions, huge charge fluctuations suddenly ap- 
pear between between weakly coupled Donor Acceptor and Catalyst, while the 
Donor-Acceptor system alone does not exhibit any fluctuations. However, the ab- 
sence of energy dissipation prevents the electron from falling on its ground-state, 
which is on the acceptor. The same model with damping shows that the electron 
finally falls on the acceptor while the catalyst has only taken transitively a fraction 
of the electronic charge (see . 

This ET is highly sensitive to small perturbations of the Donor-catalyst system 
which easily breaks the Targeted Transfer Era. We have shown for example that 
relatively small electric fields are sufficient for blocking ET at zero degree K □. 
These principles may be extended to many-site networks of electronic levels where 
the electron can choose a specific path very selectively. This path can be blocked 
and switched under small perturbations. Logical functions with one or few electrons 
could be built at the molecular level suggesting potential nanodevice applications 
and complex biological functions to be studied in living cells. These studies are left 
for further developments. 
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Fig. 8. Energies Ho (I), He (I) and Ha(I) versus electron density / for the Donor, Catalyst 
and Acceptor in the situation of Ultrafast ET. The energy barrier between Donor and Acceptor 
without Catalyst is plotted in gray. 




104 208 312 416 520 1.04 2.08 3.12 4.16 5.2 

time (ps) time (ps) 



Fig. 9. Electron density on Donor Acceptor and Catalyst versus time in the trimer in the absence 
of damping jd = ■jq = -ja = at two different time scales /ip = 2, \D = ~ b = 1) XC = b 
fiA = 1-5, XA — —0.75, Xad — ^AC = ^CD = 10 — 2 ■ The electron is initially on the Donor. 

5. Concluding Remarks 

We presented in this short paper basic principles for a non-adiabatic theory of ET. 
We briefly sketched new perspectives for understanding ultrafast ET and catalysis. 
A more complete description of this work with mathematical details, developments 
and applications shall be published elsewhere i. 

A precise example of application of our trimer model to real and well studied 
experimental problems concerns the photochemical reaction in the reaction cen- 
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Fig. 10. Same as fig-M but with damping 715 = 7c = 7A = 2 (left) or 10 (right) 



ter, where the electronic sites involved have been well identified. In this system 
the photons collected in an antenna of pigment molecules funnel to a specific site 
of the reaction center. Then, an electron is ejected and transferred within few ps 
over relatively long distances of about one nm (e.g. see ref.0 for a review). Further 
subsequent ETs follow in the biomachinery. Besides theoretical studies and numer- 
ical investigations of the molecules involved confirming the high efficiency of this 
molecular system, the Marcus theory is not sufficient for a global understanding 
of all experimental features □. We shall propose that, out of the apparent system 
complexity, our simple basic principles could help understand ultrafast ET as well 
as the puzzling features associated with it □. 
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